Eliminating inline positional errors for four-point resistance measurement

ABSTRACT

Calculating resistance correction factors includes contacting the arms of a four-arm probe with a test sample; selecting a first set of first and second arms and a second set of third and fourth arms; applying a first current from the first arm to the second arm of the first set; detecting a first voltage between the third and fourth arms of the second set; calculating a first resistance using the first voltage and current; selecting a third set of first and second arms including no more than one arm of the first set, and a fourth set of third and fourth arms including no more than one arm of the second set; applying a second current from the first arm to the second arm of the third set; detecting a second voltage between the third and fourth arms of the fourth set; calculating a second resistance using the second voltage and current; and calculating a correction factor using the first and second resistances.

CROSS-REFERENCE TO RELATED APPLICATION

This application is a divisional of co-pending application Ser. No. 12/090,425 filed Jul. 11, 2008, which is a national phase filing, under 35 U.S.C. §371(c), of International Application No. PCT/DK2006/000584, filed Oct. 17, 2006, the disclosure of which is incorporated herein by reference in its entirety.

FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not Applicable

TECHNICAL FIELD

The present invention relates to a system and method for calculation correction factors for reducing or eliminating positional errors in resistance measurements and to a method and system for performing resistance measurements for determining the sheet resistance of test samples. Furthermore, the present invention relates to a method for obtaining a characteristic electrical property of a test sample.

BACKGROUND OF THE INVENTION

Related systems and methods may be found in publications such as US 2004/0183554, U.S. Pat. No. 6,943,571, U.S. Pat. No. 4,703,252, U.S. Pat. No. 5,691,648, U.S. Pat. No. 6,747,445, US 2005/0151552, US 2005/0081609, U.S. Pat. No. 3,735,254, U.S. Pat. No. 3,456,186, WO 94/11745 and US 2005/0062448. Reference, is made to all of the above-mentioned US patent publications, all of which are hereby incorporated in the present specification by reference in their entirety for all purposes.

When performing resistance measurements for determining the resistance of a test sample, the surface of the test sample and the way that a probe used for performing the measurements are brought into contact with the surface of the test sample errors are induced by the fact that probe arms of the test probe are positioned so that the distances between individual probes vary from the assumed distances, and the current distribution in the surface of the test sample therefore is not as assumed. The present invention provides a method and system for eliminating in-line positional errors, especially in four-point probe resistance measurements.

Methods such as disclosed in US 2004/0183554 have several disadvantages, e.g. it requires knowledge of the nominal spacing. Also, there are no elimination deviations in the nominal spacing. Further the known methods are approximate. Therefore, there is a need for a method, which is exact.

SUMMARY OF THE INVENTION

According to the first aspect of the present invention, a method for calculating a correction factor for reduction of positional errors in resistance measurements using four-point probes is provided. The four-point probe has a body and four probe arms, each including a probe, the probe arms extending parallel from the body, the four-point probe including electrical contact points for establishing electrical contact between the probes and a test apparatus transmitting and receiving electrical signals The method may comprise the steps of:

-   -   positioning the probe arms in contact with a surface of the test         sample,     -   selecting a first set comprising a first and a second probe arm         and a second set comprising a third and a fourth probe arm,     -   applying a first current from the test apparatus through the         first probe arm of the first set to the second probe arm of the         first set, the first current propagating through the test         sample,     -   detecting a first induced voltage between the third and fourth         probe arms of the second set,     -   calculating a first four-point resistance, R_(first), being a         ratio of the first voltage and the first current,     -   selecting a third set comprising a first and a second probe arm,         wherein the third set includes no more than one of the probe         arms of the first set, and a fourth set comprising a third and a         fourth probe arm wherein the fourth set includes no more than         one of the probe arms of the second set,     -   applying a second current from the test apparatus through the         first probe arm of the third set to the second probe arm of the         third set, the second current propagating through the test         sample,     -   detecting a second induced voltage at the third and fourth probe         arms of the fourth set,     -   calculating a second four-point resistance, R_(second), being a         ratio of the second voltage and the second current,     -   calculating a correction factor, based on the first and second         four-point resistances.

Although the present invention is directed at four-point probes, the probes may include more than four probe arms each including at least one probe and measurements may eg. be performed using four of these probes. Also, the probe arms need not extend parallel from the probe body, but the probes, at the end of the probe arms, are preferably arranged on a line, meaning that the probe arms may have other geometrical configurations than linear or rectangular.

The probe and the test sample are preferably arranged in a test apparatus having two or more holders for holding the test sample and probe, respectively. The test probe and the test sample are brought into contact with each other, either by moving the test probe or the test sample or both so that these are brought into contact. The verification of contact between the test probe and the test sample may be verified, eg. optically/visually or using electrical methods for verifying the contact.

Ideally the surface of the test sample is flat and the probe arms of the test probe are brought into contact with the surface of the test sample so that the probes are positioned in-line having a pre-defined or known spacing. However, as this is not the case in the real world, the present invention provides a method for calculating an improved correction factor for reducing or eliminating the errors caused by the non-ideal positioning of the probes on the surface of the test sample.

In the four-point probe, a first set comprising two probe arms of the four-probe arms of the test probe is selected and a second group comprising two probe arms different from the probe arms of the first set is selected. A current is applied to one of the probe arms of the first set whereby the current is propagated through the surface of the test sample. A voltage is thereby induced in the test sample. This induced voltage may then be measured or determined using the first and second probe arms of the second set of probe arms.

From the applied current and the measured voltage a four-point resistance may be calculated as the ratio between the voltage and the current. The sheet resistance may also be measured.

A third set comprising a first and a second probe arm is selected and, provided the probe includes only four probe arms, the third set may include only one of the probe arms of the first set or alternatively the second set. A fourth set of probe arms is selected, wherein the probe arms of the fourth set include no more than one of the probe arms of the second set, or alternatively the first set.

A second current is then applied on the probe arms of the third set, or in the alternative on the probe arms of the fourth set. As the current propagates through the surface of the test sample a voltage is induced which may be determined or measured at the probe arms of the fourth set, alternatively third set. As before, a four-point sheet resistance may be calculated based on the ratio of the second voltage and the current.

In one embodiment, the calculation of the correction factor by the relation between two four-point resistances given by R_(first)/R_(second)=1−(ln(1−γ)/ln(γ)). A correction factor is then given by C_(second)=2π/ln(γ) and the sheet resistance by R_(sq)=C_(second) R_(second)

The factor γ may require the use of numerical methods. γ is a factor dependant on the geometrical configuration of the probes of a four-point probe as described elsewhere in the present specification.

According to the teachings of the present invention, the method according to the first aspect may further include performing a measurement for calculating resistance of said test sample using the correction factor.

After calculating the correction factor for a given configuration of the probe arms the probe arms are preferably not moved while performing the above-mentioned measurements and calculations. The sheet resistance may then be calculated using the correction factor as calculated above. Although the second four-point resistance is used as described above, the first four-point resistance may alternatively be used for calculating the sheet resistance of the test sample.

Advantageously, the probes may be repeatedly moved to a remote position and the steps of claim 1 and 2 are repeatedly performed for each position.

It may be necessary to determine the sheet resistance of the test sample at multiple positions and the test probe may then be stepwise moved to these locations where the calculation of the correction factor may be performed and the subsequent measurement or determination of the sheet resistance is formed.

A second aspect of the present invention relates to a method for performing resistance measurements using four-point probes using an improved correction factor for reduction or elimination of positional errors, said four-point probe having a body and four probe arms each including a probe, said probe arms extending parallel from said body, said four-point probe including electrical contact points for establishing electrical contact between said probes and a test apparatus transmitting and receiving electrical signals, said method comprising the steps of:

-   -   positioning said probe arms in contact with a surface of said         test sample,     -   selecting a first set comprising a first and a second probe arm         and a second set comprising a third and a fourth probe arm,     -   applying a first current from said test apparatus through said         first probe arm of said first set to said second probe arm of         said first set, said first current propagating through said test         sample,     -   detecting a first induced voltage at said third and fourth probe         arms of said second set,     -   calculating a first four-point resistance, R_(first), based on a         ratio of said first voltage and said first current,     -   selecting a third set comprising a first and a second probe arm,         wherein said third set includes no more than one of said probe         arms of said first set, and a fourth set comprising a third and         a fourth probe arm wherein said fourth set includes no more than         one of said probe arms of said second set,     -   applying a second current from said test apparatus through said         first probe arm of said third set to said second probe arm of         said third set, said second current propagating through said         test sample,     -   detecting a second induced voltage between said third and fourth         probe arms of said fourth set,     -   calculating a second four-point resistance, R_(second), being a         ratio of said second voltage and said second current,     -   calculating a correction factor based on said first and second         four-point resistances, and     -   calculating a resistance of said test sample using said         correction factor.

The method according to the second aspect may include any of the features mentioned in relation to the first aspect of the present invention.

A third aspect of the present invention relates to a test apparatus having a first and a second holder for holding a test sample and a test probe, respectively, a memory device for storing a computer implementation of the method according to the first and/or second aspect and a microprocessor for executing the computer implementation of the method according to the first and/or second aspect.

The test apparatus according to the third aspect of the present invention may include any of the features of the method mentioned in relation to the first and/or second aspect of the present invention.

In a fourth aspect the present invention relates to a method for obtaining a characteristic electrical property of a test sample by using a probe having a body and a plurality of probe arms each including a probe tip, the characteristic electrical property including a resistance property, the plurality being at least four, the probe arms extending parallel from the body, the probe including electrical contact points for establishing electrical contact between the probe tips and a test apparatus transmitting and receiving electrical signals to and from the probe tips, the method may comprise the following steps:

-   -   (a) positioning the probe tips in contact with a surface of the         test sample,     -   (b) selecting a subset of four probe tips for performing a         measurement by applying a current between a first set of two         probe tips and determining a voltage between a second set of two         probe tips not including any of the probe tips of the first set,         and performing the measurement for all combinations of sets of         probe tips,     -   (c) establishing a model value of the characteristic electrical         property,     -   (d) calculating model data for the characteristic electrical         property based on the model value,     -   provided the measurement and the calculated model data does not         converge, establishing an adjusted model value and repeating         step (d) using the adjusted model value,     -   provided the measurement and the calculated model data do         converge, extracting the characteristic electrical property from         the model data.

The test apparatus may comprise actuators or other devices for moving the test probe into a position where the tips are in electrical contact with the surface of the test sample. The test probe may be of a kind similar to that disclosed in publications such as US 2004/0056674 and EP 1 610 131, and the test apparatus may be of a kind similar to that disclosed in publications such as US 2005/0127929 and EP 1466 182. Reference is made to the above mentioned US patent publications, all of which is hereby incorporated into the present specification by reference.

As the probe comprises at least four probe tips, it is possible to choose combinations of two sets of two probe tips in several ways. The sets chosen in each combination do not contain any common tips. E.g. if a probe comprises the tips A, B, C and D, a first combination may be a set of A and B, and the second set being C and D. A second combination may then be a first set of A and C, and a second set of B and D. A third combination may then be a first set of B and C, and a second set of A and D. Using a probe having more tips results in further combination possibilities.

The measurements are then performed using the combinations of these sets. Model data is then calculated, based on a model value for the electrical property, for each of the combinations used. The model value may start out as an initial guess, and may be based on knowledge of the test sample.

The model data is then compared to the measured data. Provided the model data and the measured data are sufficiently close, then the method will have given an acceptable value of the characteristic electrical property of the test sample. It may be necessary to adjust the model value and repeat the calculation and comparison to the measured value or values. The adjustment may be either up or down or both if more than one characteristic electrical property is to be determined at once. The adjustment may be repeated until an acceptable result is achieved.

There may be predefined a threshold value for the difference in model data and measured data. Alternatively, the iterative process may determine a point where the difference is at a minimum.

In a further embodiment, the method may comprise calculating a corrected value of one of the measurements. One or more of the measurements may be corrected by a correction factor, or by a correction function.

In a presently preferred embodiment, the characteristic electrical property is top-layer resistance, bottom-layer resistance, resistance-area product, or any combination thereof. Other electrical properties may be determined as well.

In an advantageous embodiment of the present invention, the model data is calculated using the formula

${f_{CIPT}^{i}\left( {R_{T},R_{B},{RA}} \right)} = {\frac{R_{T}R_{B}}{2{\pi \left( {R_{T} + R_{B}} \right)}}{\left\{ {{\frac{R_{T}}{R_{B}}\left( {{K_{0}\left( \frac{w^{i}}{\lambda} \right)} + {K_{0}\left( \frac{z^{i}}{\lambda} \right)} - {K_{0}\left( \frac{y^{i}}{\lambda} \right)}} \right)} + {\ln\left( \frac{x^{i}y^{i}}{z^{i}w^{i}} \right)}} \right\}.}}$

The values x, y, z and w are distances between probe arms in a given configuration, as will be evident from the below description of the attached figures. R_(T) is the top-layer resistance, R_(B) is the bottom-layer resistance and RA is the resistance area product. K₀ is the modified Bessel function of the second kind, 0'th order. λ is the transition length given by:

$\lambda = \sqrt{\frac{RA}{R_{T} + R_{B}}}$

In an advantageous embodiment of the present invention, the method may further include performing a correction on one of said calculated model data. The correction may be performed on one or more of the calculated model data.

The method according to the fourth aspect of the present invention may include any of the features of any of the aspects one, two and/or three.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention is now to be described with reference to the attached schematic illustrations, in which:

FIG. 1 is a schematic illustration of positions of four probes and distances between them,

FIG. 2 is a schematic illustration of four probes positioned on the surface of a test sample,

FIGS. 2 a-2 c are schematic illustrations of three configurations of four-point probes and distances between the probes,

FIG. 3 is a schematic illustration of a logarithmic functions of α, β and γ as a function of b/a, with c=a,

FIG. 4 is a schematic illustration of a numerical simulation of error due to barrier parallel to the in-line direction,

FIG. 5 is a schematic illustration of numerical simulation of measurement errors on a rectangular sample,

FIG. 6 is a further illustration of simulation results, and

FIG. 7 is a block diagram representing a method according to the present invention.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1. schematically illustrates a four-point probe in contact with the surface. When a multi-point probe, such as a four-point probe, is freely suspended, the arms constituting the point probes are freely extending from the body of the multi-point probe and as such, in theory at least, are positioned in-line with the known distance between the end points of the probe and are brought into contact with the surface of the sample. The arms are displaced relative to each other partly due to surface errors or surface irregularities and possibly also due to the probes' positions relative to the surface of the sample. The spacing between the probes may vary due to uncertainty in production or vibrations. The difference in the relative spacing of the probes is illustrated in FIG. 1 in that four probe arms (10,12,14,16) are displaced relative to each other and where in between distances are illustrated by the arrows (18, 20, 22, 24,26).

Assuming that the four-point probes (10,12,14,16) are positioned in-line and with known spacing on the surface of the test sample will result in an imprecise or erroneous measurement as the current distribution and current and all voltage measured are different from the theoretical assumptions.

For an infinite sheet the four-point voltage to current ratio may be expressed as

${R_{App} = {\frac{V}{I} = {\frac{R_{B}}{2\pi}{\ln \left( \frac{xy}{wz} \right)}}}},$

where w, x, y and z are distances as given in FIG. 2 and R_(s) is the sheet resistance of the sheet. A total of 24 perturbations of interchanging electrodes exist; however these 24 perturbations can be divided into 3 groups each consisting of 8 equivalent configurations. Within each group the measured four-point resistance will differ at most by a change of sign. The four-point resistance for each group is designated R_(A), R_(B) and R_(C). These ratios are so defined that they always are positive. The following relation exists between the four-point resistances.

R _(A) =R _(B) +R _(C)

For linear configurations, representative configurations for each of these groups are illustrated in FIGS. 2 a to 2 c. The voltage-to-current ratio for these configuration are given as

$R_{A} = {\frac{V_{A}}{I_{A}} = {\frac{R_{B}}{2\pi}\ln \; \alpha}}$ ${R_{B} = {\frac{V_{B}}{I_{B}} = {\frac{R_{B}}{2\pi}\ln \; \beta}}},{R_{C} = {\frac{V_{C}}{I_{C}} = {\frac{R_{B}}{2\pi}\ln \; \gamma}}}$

where the parameters α,β and γ are given by

$\alpha = \frac{\left( {a + b} \right)\left( {b + c} \right)}{a\; c}$ $\beta = \frac{b\left( {a + b + c} \right)}{a\; c}$ $\gamma = {\frac{\left( {a + b} \right)\left( {b + c} \right)}{b\left( {a + b + c} \right)}.}$

With manipulation it can be found that the following relations exist between the parameters

$\alpha = {{\beta + 1} = {\frac{1}{\gamma - 1} = {\beta\gamma}}}$ $\beta = {{\alpha - 1} = {\frac{1}{\gamma - 1} = \frac{\gamma}{\alpha}}}$ $\gamma = {\frac{\alpha}{\alpha - 1} = {\frac{\beta + 1}{\beta} = \frac{\alpha}{\beta}}}$

Performing two measurements using the configuration illustrated in FIG. 2A and the configuration illustrated in FIG. 2C we can make the following relation

${\frac{R_{A}}{R_{C}} = \frac{{R_{B}/2}{\pi ln\alpha}}{{R_{B}/2}{\pi ln}\; \gamma}},$

and using the above relations, we can reduce this to

$\frac{R_{A}}{R_{C}} = {1 - {\frac{\ln \left( {\gamma - 1} \right)}{\ln \; \gamma}.}}$

Thus we can uniquely determine γ from the ratio R_(A) to R_(C). Due to the complexity of the relation γ must be found using numerical methods. This task may be performed using suitable routines from programmes for numerical solution of equations, eg. MINPACK. When γ has been determined we may calculate the sheet resistance from

$R_{S}^{A\; C} = {\frac{2\pi}{\ln \; \gamma}R_{C}}$

where the index AC indicates that the sheet resistance is based on the A and C configurations. In principle any other ratio of two of R_(A), R_(B) and R_(C) can be used to find any of the three parameters α, β or γ, and the sheet resistance can be found using the appropriate relation.

Since the spacings a, b and c are eliminated the method will find the correct sheet resistance independent of the real spacings. Thus any variations in positioning of the electrodes are also eliminated.

At least the following relations between four-point resistances and possible search parameters exist:

R _(A) /R _(B)=ln(α)/ln(α−1)   a1:

R _(A) /R _(B)=ln(β+1)/ln(β)   b1:

R _(A) /R _(B)=−ln(γ−1)/ln(γ)   c1:

R _(A) /R _(C)=ln(α)/(ln(α)−ln(α−1))   d1:

R _(A) /R _(C)=ln(β+1)/ln(1+1/β)   e1:

R _(A) /R _(C)=1−ln(1−γ)/ln(γ)   f1:

R _(B) /R _(C)=ln(α−1)/(ln(α)−ln(α−1))   g1:

R _(B) /R _(C)=ln(β)/(ln(β+1)−ln(β))   h1:

R _(B) /R _(C)=−ln(γ−1)/ln(γ)   i1:

Further each relation may in the implementation be represented by its reciprocal expression. In some configurations of the four-point measurement a change of sign will be needed.

The sheet resistance may be given by any of the following relations:

R _(sq) =R _(A) 2π/ln(α)   a2:

R _(sq) =R _(A) 2π/ln(β−1)   b2:

R _(sq) =R _(A) 2π/ln(γ/(γ−1))   c2:

R _(sq) =R _(B) 2π/ln(α−1)   d2:

R _(sq) =R _(B) 2π/ln(β)   e2:

R _(sq) =−R _(B) 2π/ln(γ−1)   f2:

R _(sq) =R _(C) 2π/ln(α/(α−1))   g2:

R _(sq) =R _(C) 2π/ln((β+1)/β)   h2:

R _(sq) =R _(C) 2π/ln(γ)   i2:

Any two of these relations, one from the first group and one from the second group, both containing the same parameter α, β or γ, will constitute a functional implementation of the method, e.g. a1 and g2, or e1 and b2. These relations may be derived from the previously equations. Further each of these relations may be modified by imposing R_(A)=R_(B)+R_(C).

Isolation of a search parameter in the second group and insertion in the first group will also constitute a valid implementation.

FIG. 3 is a schematic illustration of logarithmic functions of α, β and γ as a function of b/a, with c=a.

The method according to the present invention of eliminating errors due to error positioning in the in-line direction, does not increase the sensitivity to other errors such as electromagnetic radiation and positioning errors perpendicular to the in-line direction.

Errors due to electromagnetic radiation and cross-over from current source to voltage sensors e.g. due to variation in contact resistance, we can assume to result in stochastic variations directly on the measured voltage-to-current ratio. Thus we can describe the sensitivity for the normal four-point probe method simply as

${\frac{\partial R_{s}}{\partial R_{A}} = \frac{2\pi}{\ln \; \alpha}},$

and equivalently for the B and C configurations. For the presented method we have the following expression

$\frac{\partial R_{s}^{A\; C}}{\partial R_{A}} = {{{- \frac{2\pi}{\alpha \; \ln^{2}\alpha}}R_{A}\frac{\partial\alpha}{\partial R_{A}}} + {\frac{2\pi}{\ln \; \alpha}.}}$

We will use R_(A)/R_(C) to find α from

$\frac{R_{A}}{R_{C}} = {\frac{\ln \; \alpha}{\ln \left( \frac{\alpha}{\alpha - 1} \right)}.}$

By differentiating both sides with respect to R_(A) and performing manipulation we find

$\frac{\partial\alpha}{\partial R_{A}} = {\frac{\left( {\alpha^{2} - \alpha} \right)\ln^{2}\frac{\alpha}{\alpha - 1}}{R_{C}\left( {{\ln \; \alpha} + {\left( {\alpha - 1} \right)\ln \frac{\alpha}{\alpha - 1}}} \right)}.}$

By inserting equation 2 into 1 and assuming that we know α from the nominal electrode spacings, we find that the ratio

$f_{ACA} = \frac{{\partial R_{S}^{A\; C}}/{\partial R_{A}}}{{\partial R_{S}}/{\partial R_{A}}}$

is always between 0 and 1 meaning that the presented method always gives better or the same sensitivity compared to the simple method of using directly the nominal spacing. For symmetric probe spacings with a=b=c we find that the ratio f_(ACA) is 0.616, an improvement of approximately 40%. Unfortunately the ratio for the AC configuration sensitivity towards variations in R_(C) f_(ACC) is 1.894 indicating that the corrected method is almost two times as sensitive to variations in R_(C) as the uncorrected is towards variations in R_(A.) Overall the sensitivity towards variations in the voltage-to-current ratios is increased approximately by a factor of 2. This may appear to be a problem; however, since in real life measurements the measurements error is dominated by positioning errors, at least for microscopic four-point probes, the method is still an improvement over the non-corrected method.

Further it is our experience that the relative rather than the absolute noise is the same for R_(A) and R_(C.) Using that information we can write the noise in the corrected sheet resistance

$\begin{matrix} {\sigma_{R_{S}}^{{rel},c} = {\frac{1}{R_{S}}\sqrt{\begin{matrix} {{\left( {f_{ACC}R_{C}} \right)^{2}\left( \sigma_{R}^{rel} \right)^{2}} +} \\ {\left( {f_{ACA}R_{A}} \right)^{2}\left( \sigma_{R}^{rel} \right)^{2}} \end{matrix}}}} \\ {= {\sigma_{R}^{rel}\sqrt{\begin{matrix} {\left( {f_{ACC}\frac{\ln \frac{\alpha}{\alpha - 1}}{2\pi}} \right)^{2} +} \\ \left( {f_{ACA}\frac{\ln \; \alpha}{2\pi}} \right)^{2} \end{matrix}}}} \end{matrix}$

where σ_(R) _(S) ^(rel,c) and σ_(R) ^(rel) are the relative standard deviation in the corrected sheet resistance and the relative standard deviation in the measured voltage-to-current ratios, respectively. Again for a symmetric probe this factor σ_(R) _(S) ^(rel,c)/σ_(R) ^(rel) becomes 0.629. Thus we have an overall noise suppression of 37%.

The voltage-to-current ratios for standard linear four-point probe are to first order insensitive to positioning error in the direction perpendicular to the in-line direction. Thus since the correction method uses directly these ratios it too becomes insensitive to such positioning errors to first order.

Surprisingly, it turns out that the method according to the present invention eliminates the influence of non-conducting barriers parallel to the line on which the electrodes are positioned.

The above findings suggest that it might be possible to eliminate sample shape entirely. However, analytical studies with a barrier along the in-line direction show that this cannot be eliminated. On the other hand numerical studies show that the error is significantly decreased. FIG. 4 shows the simulated measurement error as a function of distance to the barrier. The distance is normalised with respect to the electrode spacings. The probe is assumed to be symmetric with a=b=c. The dashed curve shows the measurement error for a non-corrected sheet resistance while the solid curve shows the error for the corrected sheet resistance—As can be seen from the figure the method greatly improves the measurement error. For the corrected sheet resistance the error is less than 1% at a distance of just 3 times the electrode spacing. This situation occurs for the uncorrected sheet resistance at a distance of 7 times the electrode spacing.

This effect too can be eliminated if the probe is placed in the centre between two parallel barriers, or on a rectangular sample. FIG. 5 shows the error when measuring a rectangular sample.

Again we have assumed a symmetric probe. For the uncorrected sheet resistance we have used a correction factor assuming a known sample dimension and a probe located at the centre of the sample. The sample dimensions are 5 by 5 times the electrode spacing. As seen from the figure the method greatly improves the error due to misalignment of probe position. Further, unlike for the uncorrected case, it is not needed to know the exact sample dimensions if just the probe is positioned approximately in the middle of the sample.

If the sample has well-known dimensions and is very narrow, the measured value depends only weakly on the position, and the method does not provide any improvement. However, in such case the width of the sample must be known. If width is not known, the correction method is still usable to find approximate correction factors. For large samples the measurement is improved even more.

FIG. 5 is a schematic illustration of numerical simulation of measurement errors on a rectangular sample. The dashed line shows the uncorrected sheet resistance and the solid line illustrates the corrected sheet resistance.

The sample used for the simulation is a rectangular sample having a side length in the direction parallel to the in-line direction of the 5 times the electrode spacing and a side length perpendicular to the in-line direction of 5 times electrode spacing. The probe is symmetric with a=b=c. For the uncorrected sheet resistance the sample dimensions is assumed to be known and the error is calculated relative to a situation where the probe is located on the centre line.

In the following below, R_(B) denotes the bottom resistance and not as above resistance for the B-configuration for the four point measurements.

In the general case we want to find a number of parameters that characterize the sample under investigation. If we have a model that describes the expected measurement value as a function of these parameters we can search for the parameters that give the best fit to the measured data. Mathematically this is described as

${\min\limits_{P}{{{f_{Model}(P)} - f_{Measured}}}},$

where P is a set of parameter we search for, f_(Model) is the function to calculate the expected measurement values from the given model, and f_(Measured) is a vector containing the measured values. The bars represent a norm; usually the I₁ or I₂norm is chosen.

In order to decrease the influence of measurement artifacts in the parameters obtained we can introduce a correction scheme. We represent the correction scheme by the function G which takes as arguments either the measured values or the values calculated from the model. We can then write the problem of finding the model parameters as:

$\min\limits_{P}{{{{G\left\{ {f_{Model}(P)} \right\}} - {G\left\{ f_{Measured} \right\}}}}.}$

The optimal choice of G depends on the actual sample geometry and the model chosen to represent the sample. The optimal G is the one that minimizes the sensitivity for the above minimum with respect to possible measurement errors.

In the case of four-point measurement on current in plane tunnelling, CIPT, samples we have the following model

${f_{CIPT}^{i}\left( {R_{T},R_{B},{RA}} \right)} = {\frac{R_{T}R_{B}}{2{\pi \left( {R_{T} + R_{B}} \right)}}\begin{Bmatrix} {\frac{R_{T}}{R_{B}}\left( {{K_{0}\left( \frac{w^{i}}{\lambda} \right)} + {K_{0}\left( \frac{z^{i}}{\lambda} \right)} -} \right.} \\ {\left. {K_{0}\left( \frac{y^{i}}{\lambda} \right)} \right) + {\ln\left( \frac{x^{i}y^{i}}{z^{i}w^{i}} \right)}} \end{Bmatrix}}$

with the transition length

${\lambda = \sqrt{\frac{RA}{R_{T} + R_{B}}}},$

and x^(i), y^(i), z^(i) and w^(i) are the distances given in FIG. 1, for the i-th configuration. The set of parameters comprises top-layer resistance R_(T), bottom-layer resistance R_(B), and the resistance-area product RA. K₀ is the modified Bessel function of second kind, zero-th order.

We choose several groups, or sets, of contact points, each group comprising four points located approximately along a line. For each group we obtain a first four-point resistance. For a sub-set of the groups (possibly all) we obtain a second four-point resistance and calculate a corrected resistance using the method described earlier, and so on. Using the model above we calculate the expected first and second four-pointed

resistances for a first guess of R_(T), R_(B,) and RA.

We then adjust the parameters R_(T), R_(B,) and RA until we have minimized the difference from the measured and corrected data to the calculated data.

Other correction schemes may be used, which might be better for the specific sample geometry or which demands lower computational effort. The important issue is that the correction scheme is performed on both the measured data and the data calculated from the model in order to obtain comparable values.

If we imagine using a probe with 6 tips numbered 1 to 6, we can exemplify the usual measurement scheme as follows

Step 1: Measure the four-point resistance using

f _(Measured) ¹ =R _(4pp)(n _(i1)=1,n _(v1)=2,n _(v2)=3,n _(i2)=4)

f _(Measured) ² =R _(4pp)(n _(i1)=2,n _(v1)=3,n _(v2)=4,n _(i2)=5)

f _(Measured) ³ =R _(4pp)(n _(i1)=3,n _(v1)=4,n _(v2)=5,n _(i2)=6)

f _(Measured) ⁴ =R _(4pp)(n _(i1)=1,n _(v1)=3,n _(v2)=4,n _(i2)=5)

f _(Measured) ⁵ =R _(4pp)(n _(i1)=1,n _(v1)=3,n _(v2)=5,n _(i2)=6)

where n_(i1) is the tip number where current is applied, n_(i2) is the tip number where current drawn out, n_(v1) and n_(v1) are the tips over which the voltage is measured.

Step 2: Guess initial parameters R_(T)=R_(T,0), R_(B)=R_(B,0), RA=RA₀.

Step 3: Calculate model data using the equations:

f _(Model) ¹ =f _(CIPT)(x ¹ =d ₁₃ ,y ¹ =d ₂₄ ,z ¹ =d ₁₂ ,w ¹ =d ₃₄)

f _(Model) ² =f _(CIPT)(x ² =d ₂₄ ,y ² =d ₃₅ ,z ² =d ₂₃ ,w ² =d ₄₅)

f _(Model) ³ =f _(CIPT)(x ³ =d ₃₅ ,y ³ =d ₄₅ ,z ³ =d ₃₄ ,w ³ =d ₅₆)

f _(Model) ⁴ =f _(CIPT)(x ⁴ =d ₁₄ ,y ⁴ =d ₃₅ ,z ⁴ =d ₁₃ ,w ⁴ =d ₄₅)

f _(Model) ⁵ =f _(CIPT)(x ⁵ =d ₁₅ ,y ⁵ =d ₃₆ ,z ⁵ =d ₁₃ ,w ⁵ =d ₅₆)

where d_(kl) is the distance between the k-th and the l-th tip.

Step 4: Calculate the objective function

F _(Obj)=Σ_(i=1) ⁵(f _(Model) ^(i) −f _(Measured) ^(a))²

where we have used the I₂-norm in the equation

$\min\limits_{P}{{{{G\left\{ {f_{Model}(P)} \right\}} - {G\left\{ f_{Measured} \right\}}}}.}$

Step 5: Adjust parameters and redo steps 3 and 4 until F_(obj) have reached minimum.

Step 6: Extract resulting parameters R_(T), R_(B) and RA.

The following is an example where correction is performed on two sets of measurements. In the case of using correction on the first two sets, or groups, of tips the scheme becomes

Step 1: Measure the four-point resistance using

f _(Measured) ^(1A) =R _(4pp)(n _(i1)=1,n _(v1)=2,n _(v2)=3,n _(i2)=4)

f _(Measured) ^(1C) =R _(4pp)(n _(i1)=1,n _(v1)=4,n _(v2)=3,n _(i2)=2)

f _(Measured) ^(2A) =R _(4pp)(n _(i1)=2,n _(v1)=3,n _(v2)=4,n _(i2)=5)

f _(Measured) ^(2C) =R _(4pp)(n _(i1)=2,n _(v1)=5,n _(v2)=4,n _(i2)=3)

f _(Measured) ³ =R _(4pp)(n _(i1)=3,n _(v1)=4,n _(v2)=5,n _(i2)=6)

f _(Measured) ⁴ =R _(4pp)(n _(i1)=1,n _(v1)=3,n _(v2)=4,n _(i2)=5)

f _(Measured) ⁵ =R _(4pp)(n _(i1)=1,n _(v1)=3,n _(v2)=5,n _(i2)=6)

Step 2: perform corrections:

G _(Measured) ¹ =G _(Sheet)(f _(Measured) ^(1A) ,f _(Measured) ^(1C))

G _(Measured) ² =G _(Sheet)(f _(Measured) ^(2A) ,f _(Measured) ^(2C))

G_(Measured) ³=f_(Measured) ³

G_(Measured) ⁴=f_(Measured) ⁴

G_(Measured) ⁵=f_(Measured) ⁵

Where G_(sheet) represents the method described previously.

Step 3: Guess the initial parameters

R_(T)=R_(T,0), R_(B)=R_(B,0) RA=RA₀

Step 4: Calculate model data using:

f _(Model) ^(1A)(x ^(1A) =d ₁₃ ,y ^(1A) =d ₂₄ ,z ^(1A) =d ₁₂ ,w ^(1A) =d ₃₄)

f _(Model) ^(1C)(x ^(1C) =d ₁₃ ,y ^(1C) =d ₂₄ ,z ^(1C) =d ₂₃ ,w ^(1C) =d ₁₄)

f _(Model) ^(2A)(x ^(2A) =d ₂₄ ,y ^(2A) =d ₃₅ ,z ^(2A) =d ₂₃ ,w ^(2A) =d ₄₅)

f _(Model) ^(2C)(x ^(2C) =d ₂₄ ,y ^(2C) =d ₃₅ ,z ^(2C) =d ₃₄ ,w ^(2C) =d ₂₅)

f _(Model) ³(x ³ =d ₃₅ ,y ³ =d ₄₅ ,z ³ =d ₃₄ ,w ³ =d ₅₆)

f _(Model) ⁴(x ⁴ =d ₁₄ ,y ⁴ =d ₃₅ ,z ⁴ =d ₁₃ ,w ⁴ =d ₄₅)

f _(Model) ⁵(x ⁵ =d ₁₅ ,y ⁵ =d ₃₅ ,z ⁵ =d ₁₃ ,w ⁵ =d ₅₆)

Step 5: Perform pseudo correction on calculated data:

G _(Model) ¹ =G _(Sheet)(f _(Model) ^(1A) ,f _(Model) ^(1C))

G _(Model) ² =G _(Sheet)(f _(Model) ^(2A) ,f _(Model) ^(2C))

G_(Model) ³=f_(Model) ³

G_(Model) ⁴=f_(Model) ⁴

G_(Model) ⁵=f_(Model) ⁵

Step 6: Calculate the objective function

F _(Obj)=Σ_(i=1) ⁵(G _(Model) ^(i) −G _(Measured) ^(i))²

Where the I₂ norm is used in the equation

$\min\limits_{P}{{{{G\left\{ {f_{Model}(P)} \right\}} - {G\left\{ f_{Measured} \right\}}}}.}$

Step 7: Adjust parameters and redo steps 4, 5 and 6 until F_(obj) have reached minimum.

Step 8: Extract resulting parameters R_(T), R_(B) and RA.

The method is illustrated in FIG. 6

FIG. 7 illustrates a simulation of the performance of the correction scheme, f_(Measured), as been simulated using the configuration of a real 12 point probe with smallest tip-spacing of 1 μm. For each simulation 8 sets or groups of tips have been used. The simulation uses R_(T)=1Ω, R_(B)=1Ω and RA=5 Ω/μm². Each tip position varies with a random square distribution of ±50 nm. The figure clearly shows that the variation in the extracted RA decreases with increasing level of correction at least until the set with the 3rd smallest spacing.

TABLE 1 No. of corrections 0 1 2 3 Relative standard deviation R_(T) (%) 700 120 4.1 4.8 Relative standard deviation R_(B) (%) 51.2 26.7 3.9 4.2 Relative standard deviation RA (%) 685 61.3 10.2 12.5 Mean R_(T) 2.7 * 10⁶ 2.02 1.00 1.00 Mean R_(B) 1.05 0.98 1.00 1.00 Mean RA 900 1.32 1.01 1.01 R_(T) = 1 Ω, R_(B) = 1 Ω, RA = 1 Ω/μm²

TABLE 2 No. of corrections 0 1 2 3 Relative standard deviation R_(T) (%) 1.0 * 10³ 29.4 2.0 1.8 Relative standard deviation R_(B) (%) 29.2 9.8 1.8 1.1 Relative standard deviation RA (%) 1.0 * 10³ 37.6 7.3 6.8 Mean R_(T) 1.1 * 10⁷ 1.06 1.00 1.00 Mean R_(B)  0.96 0.99 1.00 1.00 Mean RA 2.4 * 10⁵ 2.14 1.99 2.02 R_(T) = 1 Ω, R_(B) = 1 Ω, RA = 2 Ω/μm²

TABLE 3 No. of corrections 0 1 2 3 Relative standard deviation R_(T) (%) 9.9 2.4 0.8 0.5 Relative standard deviation R_(B) (%) 7.9 1.5 0.9 0.9 Relative standard deviation RA (%) 27.2 18.7 5.1 2.7 Mean R_(T) 1.02 1.00 1.00 1.00 Mean R_(B) 1.00 1.00 1.00 1.00 Mean RA 5.14 5.20 5.05 4.98 R_(T) = 1 Ω, R_(B) = 1 Ω, RA = 5 Ω/μm²

TABLE 4 No. of corrections 0 1 2 3 Relative standard deviation R_(T) (%) 5.8 0.8 0.4 0.2 Relative standard deviation R_(B) (%) 4.7 1.1 0.8 0.9 Relative standard deviation RA (%) 21.1 13.3 3.9 1.5 Mean R_(T) 0.99 1.00 1.00 1.00 Mean R_(B) 1.01 1.00 1.00 1.00 Mean RA 10.6 10.3 10.0 10.0 R_(T) = 1 Ω, R_(B) = 1 Ω, RA = 10 Ω/μm²

TABLE 5 No. of corrections 0 1 2 3 Relative standard deviation R_(T) (%) 3.8 0.1 0.1 0.1 Relative standard deviation R_(B) (%) 18.8 10.0 5.9 2.0 Relative standard deviation 34.5 11.2 7.7 2.2 RA (%) Mean R_(T) 1.00 1.00 1.00 1.00 Mean R_(B) 0.96 0.97 0.99 1.00 Mean RA 110 103 101 100 R_(T) = 1 Ω, R_(B) = 1 Ω, RA = 100 Ω/μm²

Tables 1 to 5 summarise simulations using different values of RA. The results show that the correction method generally decreases the error on the extracted parameters and lower the value of RA where meaningful parameters can be extracted.

The magneto resistance of the CIPT-structure is not included in the demonstration of the method, since it does not change the basic concept. The method can be applied both with subsets of a multipoint probe and measurements using several multipoint probes. The correction can be applied on any number of the sets of tips.

The correction need not to be applied on the smallest spacing, but can be applied on any of the sets or groups of tips. The method can be applied on other geometries than the CIPT structure provided that the model is changed accordingly. Other correction schemes than the presented may be used. The crucial point is that correction is applied on both the measured and calculated values, to create comparable data. Calculation of corrected values may be an integrated part of performing the measurement. New measurements may be performed after an initial correction in order to refine the result. As described previously any of the linear independent configurations may be used in the correction. And the search parameter, i.e. model value, may be any of α, β or γ as described earlier.

FIG. 6 schematically illustrates variations in fitted RA with real RA=5 Ω/μm². In the figure are shown results with no correction 27 a, correction on smallest spacing 27 b, correction on 2 smallest spacings 27 c and correction on 3 smallest spacings 27 d.

FIG. 7 schematically illustrates steps of an embodiment of the method according to the teachings of the present invention. The starting point is step 28. Sets of four-point measurement are performed in step 30. Prior to the measurement being performed a probe has been positioned in electrical contact with a surface of a test sample.

For a subset of the measurements, alternative linear independent configurations are measured in step 32. The measurements may be performed in different combinations, i.e. there are multiple possibilities where the current may be applied, and likewise for the choice of points where the corresponding voltage is to be measured.

Calculate corrected resistance for the subset in step 34. After performing the measurements in steps 30 and 32, correction values are then calculated.

Choose a first guess for R_(T), R_(B) and RA in step 36. The first guess for the electrical properties to be determined is established.

Calculate expected resistances for all configurations in step 38. Based on the first guess for R_(T), R_(B) and RA the expected, or theoretical, model value is calculated.

Perform correction for the subset of configurations in step 40. The model values are corrected in the same way as the measured values.

Compare measured corrected resistances to calculated corrected resistances in step 42. A comparison is made between the corrected, measured values and the corrected, model values.

In step 44 a check is performed to see if the corrected, measured values and the corrected, model values converge. If the corrected, measured values and the corrected, model values are sufficiently close, the method has yielded an acceptable answer.

Output R_(T), R_(B) and RA 46, and the method ends in step 48.

The values may e.g. be outputted to a recording device or a printer device.

If the corrected, measured values and the corrected, model values do not converge, the guess for R_(T), R_(B) and RA is adjusted in step 50. The adjusted values are then returned for a new calculation in steps 38-44. The process may be iterated until a sufficiently acceptable answer is obtained. The acceptance of the values R_(T), R_(B) and RA may be based on determining a minimum between the corrected, measured values and the corrected, model values.

The present invention may be characterised by the following points:

1. A method for calculating a correction factor for reduction of positional errors in resistance measurements using four-point probes, said four-point probe having a body and four probe arms each including a probe, said probe arms extending parallel from said body, said four-point probe including electrical contact points for establishing electrical contact between said probes and a test apparatus transmitting and receiving electrical signals, said method comprising the steps of:

-   -   positioning said probe arms in contact with a surface of said         test sample,     -   selecting a first set comprising a first and a second probe arm         and a second set comprising a third and a fourth probe arm,     -   applying a first current from said test apparatus through said         first probe arm of said first set to said second probe arm of         said first set, said first current propagating through said test         sample,     -   detecting a first induced voltage between said third and fourth         probe arms of said second set,     -   calculating a first four-point resistance, R_(first), being a         ratio of said first voltage and said first current,     -   selecting a third set comprising a first and a second probe arm,         wherein said third set includes no more than one of said probe         arms of said first set, and a fourth set comprising a third and         a fourth probe arm wherein said fourth set includes no more than         one of said probe arms of said second set;     -   applying a second current from said test apparatus through said         first probe arm of said third set to said second probe arm of         said third set, said second current propagating through said         test sample,     -   detecting a second induced voltage at said third and fourth         probe arms of said fourth set,     -   calculating a second four-point resistance, R_(second), being a         ratio of said second voltage and said second current,     -   calculating a correction factor based on said first and second         four-point resistance.

2. The method according to point 1, wherein said method further includes:

-   -   performing a measurement for calculating resistance of said test         sample using said correction factor.

3. The method according to point 2, wherein said probe is repeatedly moved to a remote position, and said steps of points 1 and 2 are repeatedly preformed for each position.

4. A method for performing resistance measurements using four-point probes using a correction factor for reduction or elimination of positional errors, said four-point probe having a body and four probe arms each including a probe, said probe arms extending parallel from said body, said four-point probe including electrical contact points for establishing electrical contact between said probes and a test apparatus transmitting and receiving electrical signals, said method comprising the steps of:

-   -   positioning said probe arms in contact with a surface of said         test sample,     -   selecting a first set comprising a first and a second probe arm         and a second set comprising a third and a fourth probe arm,     -   applying a first current from said test apparatus through said         first probe arm of said first set to said second probe arm of         said first set, said first current propagating through said test         sample,     -   detecting a first induced voltage at said third and fourth probe         arms of said second set,     -   calculating a first four-point resistance, R_(first), being a         ratio of said first voltage and said first current,     -   selecting a third set comprising a first and a second probe arm,         wherein said third set includes no more than one of said probe         arms of said first set, and a fourth set comprising a third and         a fourth probe arm wherein said fourth set includes no more than         one of said probe arms of said second set,     -   applying a second current from said test apparatus through said         first probe arm of said third set to said second probe arm of         said third set, said second current propagating through said         test sample,     -   detecting a second induced voltage between the third and fourth         probe arms of said fourth set.     -   calculating a second four-point resistance, R_(second), being a         ratio of said second voltage and said second current,     -   calculating a correction factor based on said first and second         four-point resistances     -   calculating resistance of said test sample using said correction         factor.

5. The method according to point 4, wherein said probe is repeatedly moved to a remote position, and said steps of point 4 are repeatedly performed for each position.

6. A test apparatus having a first and a second holder for holding a test sample and a test probe, respectively, a memory device for storing a computer implementation of the method according to point 1 and/or point 4 and a microprocessor for executing the computer implementation of the method according to point 1 and/or point 4. 

1-6. (canceled)
 7. A method for obtaining a characteristic electrical property of a test sample by using a probe having a body and at least four probe arms each including a probe tip, said characteristic electrical property including a resistance property, said probe arms extending parallel from said body, said probe including electrical contact points for establishing electrical contact between said probe tips and a test apparatus operable for transmitting and receiving electrical signals to and from said probe tips, said method comprising the following steps: (a) positioning said probe tips in contact with a surface of a test sample, (b) selecting a subset of four probe tips for performing a measurement by applying a current between a first set of two probe tips and determining a voltage between a second set of two probe tips not including any of the probe tips of said first set, and performing said measurement for all combinations of sets of probe tips, (c) establishing a model value of said characteristic electrical property, (d) calculating model data for said characteristic electrical property based on said model value, (e) if said measurement and said calculated model data do not converge, establishing an adjusted model value and repeating step (d) using said adjusted model value, and (f) if said measurement and said calculated model data do converge, extracting said characteristic electrical property from said model data.
 8. The method according to claim 7, wherein said step (b) further comprises calculating a corrected value of one of said measurements.
 9. The method according to claim 7, wherein said characteristic electrical property is selected from the group consisting of top-layer resistance, bottom-layer resistance, resistance-area product, and any combination thereof.
 10. The method according to claim 7, wherein said model data is calculated using the formula ${f_{CIPT}^{i}\left( {R_{T},R_{B},{RA}} \right)} = {\frac{R_{T}R_{B}}{2{\pi \left( {R_{T} + R_{B}} \right)}}{\begin{Bmatrix} {\frac{R_{T}}{R_{B}}\left( {{K_{0}\left( \frac{w^{i}}{\lambda} \right)} + {K_{0}\left( \frac{z^{i}}{\lambda} \right)} -} \right.} \\ {\left. {K_{0}\left( \frac{y^{i}}{\lambda} \right)} \right) + {\ln\left( \frac{x^{i}y^{i}}{z^{i}w^{i}} \right)}} \end{Bmatrix}.}}$
 11. The method according to claim 7, wherein said method further includes performing a correction on one of said calculated model data.
 12. (canceled) 